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Ultraviolet Fen emission in z ~ 2 quasars 
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ABSTRACT 

We present spectra of six luminous quasars at z ~ 2, covering rest wavelengths 
1600-3200 A. The fluxes of the UV Fen emission lines and Mgn A2798 doublet, the 
line widths of Mgn, and the 3000 A luminosity were obtained from the spectra. These 
quantities were compared with those of low-redshift quasars at z = 0.06 — 0.55 studied 
by Tsuzuki et al. In a plot of the Fen(UV)/Mgn flux ratio as a function of the cental 
black hole mass, Fen(UV)/Mgn in our z ~ 2 quasars is systematically greater than 
in the low-redshift quasars. We confermed that luminosity is not responsible for this 
excess. It is unclear whether this excess is caused by rich Fe abundance at z ~ 2 over 
low-redshift or by non-abundance effects such as high gas density, strong radiation 
field, and high microturbulent velocity. 

Key words: galaxies: abundances - galaxies: active - line: formation - quasars: 
emission lines 



1 INTRODUCTION 

According to the models of explosive nucleosynthesis, much 
of the iron comes from Type la supernovae, while a elements 
such as O and Mg come from Type II supernovae. Because 
the difference in lifetime of the progenitors, it is generally 
considered that the iron enrichment delays relative to a ele- 
ments by 1-2 billion years (Hamann & Ferland 1993; Yoshii 
et al. 1996, 1998). If Fell/Mgll, the relative strengths of 
Fell emission lines and the Mgll A2798 doublet, reflects the 
Fe/Mg abundance ratio, there will be a break in Fell/Mgll 
at high redshift. Despite of much efforts made by many ob- 
servational groups (e.g., Elston et al. 1994; Kawara et al. 
1996; Dietrich et al. 2002 , 2003; Iwamuro et al. 2002, 2004; 
Freudling et al. 2003; Maiolino et al. 2003; Tsuzuki et al. 
2006; Matsuoka et al. 2007, 2008a; Kurk et al. 2007), there 
have been found no signs of such a break; Fell/Mgll looks 
constant from low-redshift up to z ~ 6.5 with large scatter. 

No break in Fell/Mgll might reflect a significantly 
shorter delay-time of 0.2—0.6 Gyr, as suggested by Friaga & 
Terlevich (1998), Matteucci & Recchi (2001), and Granato 
et al. (2004). The expected break can also be obscured by 
non-abundance effects. Simulations of Fe II emitting regions, 
assuming either photoionization or shocks, imply that the 
Fe abundance is not only one parameter which controls the 
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Fell strength, but several non-abundance factors can also af- 
fect it. Such non-abundance factors include spectral energy 
distribution (SED) of the central source, strength of the ra- 
diation field, and the gas density of Broad Emission Line 
Region (BELR) clouds. Recently, Verner et al. (2003) and 
Baldwin et al. (2004) pointed out that a large microturbu- 
lence velocity may be responsible for strong Fell emission. 
Tsuzuki et al. (2006) have studied non-abundance factors 
by using spectra of a low-redshift sample of 14 quasars, cov- 
ering wide rest-wavlengths 1000—7300 A, and claimed that 
the Fell strength correlates with the mass of the central 
black hole, the line width, and the X-ray photon index. 

In this paper, we present spectra of six quasars at 
z ~ 2, and compare with those in the low-redshift sample. 
Throughout this paper, a cosmology with £l m = 0.3, Qa = 
0.7, and Ho = 70 km s _1 Mpc" 1 is assumed. 



2 OBSERVATIONS 

Six quasars were selected for optical spectroscopy from the 
catalog by Veron-Cetty & Veron (2003). According to the 

catalog, these are luminous with Mb = —28 31 at z — 

2.0 - 2.3 , bright enough to take optical Fell emission lines 
through near-infrared spectroscopy at a later opportunity. 

GMOS on Gemini-South Telescope are used in the long- 
slit mode with grating R150_G5326 and order sorting filter 
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Figure 1. Spectra of six quasars. The shaded area indicates the 
region where the spectra are affected by the telluric absorption. 

OG515_G0330. Wavelengths observed are in 5400-9800 A, 
corresponding to rest wavelengths 1800—3300 A at z ~ 2 
quasars. The slit width is TO" and the spectral resolution is 
3.286 A pixel" 1 . The grating was centered at 8150 A for the 
first three exposures and changed to 8250 A for the following 
exposures, filling up the gaps of the CCD chip array. Wave- 
lengths were calibrated using the CuAr arc lamp taken at the 
both central wavelengths of the grating. LTT1788 (Hamuy 
et al. 1994) was used for flux-scaling. The observing log is 
summarlized in Table 1. 

Individual spectral frames were processed using the 
Gemini IRAF package verl.9.1. Sky of individual frames 
is subtracted with GSSKYSUB, and then SCOMBINE was 
used to combine into the final spectral frame. The correc- 
tions for telluric absorption were not applied. Reduced spec- 
tra are shown in Figure 1 . Redshifts are determined from fit 
of Mgll emission line. With these measured redshifts, ab- 
solute B magnitudes are calculated assuming a cosmology 
with D. m = 0.3, ft A = 0.7, and H = 70 km s _1 Mpc -1 . 
Redshifts and absolute B magnitudes are also listed in Ta- 
ble 1. 



3 MEASUREMENT OF EMISSION LINES 

Prior to measuring physical quantities such as Fe II emission 
lines, the quasar spectra were dereddened for the Galac- 
tic extinction according to the dust map by Schlegel et al. 
(1998) using the Milky Way extinction curve by Pei et al. 
(1992). Eb-v of the Galactic extinction is listed in Table 1. 



In the shaded area in Figure 1, the telluric absorption fea- 
tures are seen. We have not applied any correction for the 
telluric absorption. Instead, the intensities in the shaded ar- 
eas were estimated by fitting a linear function to assumed 
data points locating on either side free from the telluric fea- 
tures. 

3.1 Fell UV emission lines 

Fell emission lines are heavily blended with each other, 
forming the broad features from 2000—3000 A. It is desirable 
to observe a wide range of wavelengths in such a way that 
the power-law and Balmer continua are accurately deter- 
mined as made by Tsuzuki et al. (2006) . However, observing 
such a wide range is not feasible in most cases. In fact, the 
present observations are limited to a rest wavelength range 
from 1600-3200 A. 

We applied a simple alternative in which a linear func- 
tion is fit to the data in rest wavelengths 2190—2230 A and 
2660—2700 A. Differences between the spectrum and the re- 
sultant best-fit function, which are marked as shade in Fig- 
ure 2, are summed up in a wavelengh range of 2240-2650 A. 
The summed- up differences, as denoted by /(2240— 2650 A), 
should contribute significant part of Fe//(2000 - 3000 A) 
which is the total Fe II emission line flux in 2000—3000 A. To 
check the relationship between them, we have applied this 
alternative to the low-redshift quasars studied by Tsuzuki 
et al. (2006) for which Fez/ (2000 - 3000 A) is known. The 
results are shown in Figure 3(a). This figure indicates that 
the relation is linear and /(2240 — 2650 A) is approximately 
40% of Fe 7/(2000 - 3000 A). A least-squares fitting to the 
data gives the following relation: 

logFe//=log/ + 0.402(±0.142) (1) 

Here Fe 11 = Fe II (2000 - 3000 A) and / = /(2240-2650 A). 

This relation will be used to convert observed /(2240 — 
2650 A) to the total flux of the Fell emission lines in 
2000—3000 A in the later part of this paper. 

3.2 Mgll emssion lines 

To measure the flux and the full width at half maximum 
(FWHM) of the Mgll A2798 doublet, Tsuzuki et al. (2006) 
fitted a single Gaussian component to the spectrum where 
the power- low and Balmer continua, and Fe II emission fea- 
tures were already subtracted. Again, we are not allowed to 
apply their method because of our limited wavelength range. 

Our alternative is illustrated in Figure 4. A linear func- 
tion is fit to the data in 2660-2700 A and 2930-2970 A 
where contributions from the Fell and Mgll emission lines 
are relatively weak and thus the power-law continuum can 
be defined 1 . This fitted function is subtracted from the spec- 
trum, as shown in Figure 4. To measure the Mgll FWHM, 
we first applied smoothing to that subtracted spectrum, 
then measured the velocity range within which the flux be- 
come more than half of its maximum value and defined it 
as FWHM. To minimize contributions from Fell emission, 
we only integrate the flux within -FWHM < u(Mg II) < 

1 3000-3050 A would be a better choice than 2930-2970 A if the 
CCD fringe pattern can be well removed in these wavelenghts. 
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Table 1. Observing Log for Gemini Quasars 



Object 


a 


5 


Redshift a 


M B b 


Eb-v c 


Exposure Time 


Date 










[mag] 


[mag] 


[s] 




B0226-104 


02 28 39.2 


-10 11 10 


2.276 


-29.7 


0.03 


120 


2004 Sep 18 


B0421+019 


04 24 08.6 


+02 04 25 


2.059 


-27.8 


0.19 


600 


2004 Sep 18 


CTQ254 


04 30 14.6 


-36 26 47 


2.118 


-27.7 


0.02 


1500 


2004 Sep 18 


FIRSTJ2149-0811 


21 49 48.2 


-08 11 16 


2.128 


-28.9 


0.04 


150 


2004 Sep 20 


LBQS2209-1842 


22 12 10.4 


-18 27 38 


2.093 


-27.4 


0.03 


960 


2004 Sep 20 


PHL424 


23 13 24.5 


+00 34 45 


2.087 


-28.5 


0.04 


300 


2004 Sep 21 



a Redshift as determined from fit of Mg II emission line. 

k Absolute B magnitude with f2 m = 0.3, Qa = 0.7, and Ho = 70 km s _1 Mpc -1 . 
c Galactic extinction Eg_y taken from Schlegel et al. (1998) 
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Figure 2. Measuring the UV Fe II line flux. A linear fuction is fit 
to the data in the continuum windows (thick bars). The best-fit 
function is indicated by the dashed line. The shaded area indicate 
/(2240-2650 A), which is then converted to Fe 7/(2000 -3000 A) 
using equation (1). 
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Figure 3. Comparison of the physical quantities in the low- 
redshift sample measured between the two methods; data along 
the horizontal- axis are measured by our method and those along 
the vertiacal-a,xis by Tsuzuki ct al. (2006). The dashed line in- 
dicate y = x where two measuements result in the identical 
values, and the solid line indicates the best-fit line. Each panel 
shows (a) Fe 11(2000-3000 A) [ergs s" 1 cm" 2 ], (b) Mgll(total) 
[ergs s- 1 cm" 2 ], (c) FWHM(Mgll) [km s" 1 ], (d) 3000 A Lumi- 
nosity [10 37 W]. 



FWHM and define it as /(Mgll < FWHM). Again, we used 
the low-redshift sample by Tsuzuki et al. (2006) to check 
that the relation between /(Mgll < FWHM) and the to- 
tal Mgll flux Mg Il(total). The results are shown in Figure 
3(6) - (c). /(Mgll < FWHM) has linear relation to the total 
Mgll line flux measured by Tsuzuki et al. (2006), and the 
least-squares best-fit to the data gives the following relation: 

log Mgii(total)= log /(Mgll < FWHM) - 0.049(±0.084)(2) 

This will be used to obtain Mg II (total) in the later part 
of this paper. It is noted that there are no significant dif- 
ferences in FWHM of Mgll between the measurements by 
Tsuzuki et al. (2006) and our alternative. 



3.3 Luminosity 

McLure & Jarvis (2002) gives a method for estimating black 
hole masses of quasars using the FWHM of the Mg II emis- 
sion line and the continuum luminosity at 3000 A. The equa- 
tion is as follows: 



M(7 



= 3.37 



/ AL3000 \ 

uo 37 wJ 



FWHM(Mg II) 



km s 



(3) 



Note that the uncertainty of equation (3) is a factor of 2.5. 

Hence the monochromatic luminosity AL3000 has to be 
measured in order to estimate the black hole mass. Unfortu- 
nately, however, our spectra are significantly affected by tel- 
luric absorption at 3000 A. We thus extrapolated the linear 
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Figure 4. Measuring the Mg II flux. Since there are Fe II emis- 
sions underneath Mgll, we measured the flux only over the ve- 
locity range — FWHM < t>(Mg II) < FWHM and defined it as 
/(Mgll < FWHM) (shaded area). 



by Tsuzuki et al. (2006). Our quasars have an absolute lu- 
minosity of Mb < —27.4, which are much luminous than 
the low-redshift quasars having an absolute luminosity of 
M B > -26.3. 

As can be seen in Figure 5(a), the five z ~ 2 quasars do 
not follow the correlation found by Tsuzuki et al. (2006). All 
of them have Fe Il(UV)/Mg II greater than expected from 
the Tsuzuki's correlation. What is the cause of the large 
FeIl(UV)/MgII value in the z ~ 2 quasars relative to the 
low-redshift quasars? Because the five z ~ 2 quasars are 
much luminous than the 14 low-redshift quasars, the lumi- 
nosity effect is examined. As shown in Figure 5(6), the lumi- 
nosity effect is not responsible for the large FeIl(UV)/MgII 
value in the z ~ 2 quasars. The real cause would be evo- 
lution in FeIl(UV)/MgII or non-abundance effects such as 
the sepectral energy distribution of the continuum from the 
central source, the strength of the radiation field and the 
gas density of BLR clouds as well as the microturbulence of 
BLR gas (Verner et al. 2003; Baldwin et al. 2004). Further 
investigations are required using large samples of quasars. 
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function, which was used to measure /(Mgll < FWHM), to 
3000 A and read the value at 3000 A as a monochromatic 
flux. Figure 3(d) compares the 3000 A luminosity AL3000 
measured by Tsuzuki et al. (2006) and our extrapolated 3000 
A luminosity A/3000 • Again, the linear relation was obtained 
as follows: 

log AL3000 = log A/3000 - 0.112(±0.079) (4) 

This will be used to convert extrapolated luminosities 
to the real luminosity in the later part of this paper. 



4 RESULTS 

In Table 2, the Fell emission line flux Fe//(2000 - 3000 A), 
the Mgll line flux Mgll (total), the FWHM (Mgll), the 3000 
A luminosity AL3000, and the blackhole mass Mbh derived 
from equation (3), are given. Note that, for B0226-104, Mgll 
emission line is heavily affected by telluric absorption line, 
so that we did not measure Mgll flux and Mgll FWHM. 
A plot of Fe 7/(2000 - 3000 A)/ Mgll (total) against Mbh is 
given in Figure 5(a) and 3000 A luminosity AL3000 in Figure 
5(6). 



5 DISCUSSION 

Analyzing 14 low-redshift quasars, Tsuzuki et al. (2006) 
found the correlation between the flux ratio FeIl(UV)/MgII 
and the blackhole mass. This relation is shown by the dotted 
line in Figure 5(a). Filled circles are our quasars at z ~ 2.0 
and open circles are low-redshift quasars at z = 0.06 — 0.55 
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Table 2. Measured Physical Quantities 



Object 


Fell (2000-3000 A) a 


Mgll(total) 


FWHM(Mg 


Il) b 


AL3000 








[10~ 14 ergs s~ 1 cur 2 ] 


[10~ 14 ergs s" 1 cm- 2 ] 


[km/s] 




[10 37 W] 


[10 9 M Q ] 


B0226-104 


19.3 (+7.5/-5.4) 














B0421+019 


5.46 (+2.1/-1.5) 


1.24 (+0.26/-0.22) 


4410 (+520/- 


-470) 


369 (+74/- 


62) 


1.05 + 0.269 


CTQ254 


3.11 (+1.2/-0.87) 


0.726 (+0.16/-0.13) 


4750 (+570/- 


-510) 


161 (+32/- 


27) 


0.828 + 0.212 


FIRSTJ2149-0811 


14.1 (+5.5/-3.9) 


2.16 (+0.46/-0.38) 


3900 (+460/- 


-410) 


500 (+100/- 


-83) 


0.951 + 0.243 


LBQS2209-1842 


8.78 (+3.4/-2.5) 


1.70 (+0.36/-0.30) 


4180 (+500/- 


-440) 


360 (+72/- 


60) 


0.935 + 0.239 


PHL424 


19.5 (+7.6/-5.5) 


2.81 (+0.60/-0.49) 


5000 (+590/- 


-530) 


587 (+120/- 


-98) 


1.69 + 0.431 



a Fell is dinned in a wavelength range of 2000-3000 A. 
b FWHM of the Mg II emission line. 
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(a) Fell/Mgll vs Black Hole Mass 
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(b) Fell/Mgll vs Luminosity 
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Figure 5. (a) Relation between Fe Il(UV)/Mg II flux ratio and black hole mass. Filled circles are our quasars at 2 ~ 2.0. Open circles 
are low-redshift quasars at z = 0.06 — 0.55 by Tsuzuki et al. (2006). Dotted line indicates the correlation found in Tsuzuki et al. (2006). 
(b) Same as in (a), but for 3000 A luminosity. 



McLure, R. J., & Jarvis, M. J. 2002, MNRAS, 337, 109 
Pei, Y. C. 1992, ApJ, 395, 130 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 
500, 525 

Tsuzuki, Y., Kawara, K., Yoshii, Y., Oyabu, S., Tanabe, 

T., & Matsuoka, Y. 2006, ApJ, 650, 57 
Verner, E., Bruhweiler, F., Verner, D., Johansson, S., & 

Gull, T. 2003, ApJ, 592, L59 
Veron-Cetty, M.-R, & Veron, P. 2003, A&A, 412, 399 
Yoshii, Y., Tsujimoto, T., & Kawara, K. 1998, ApJ, 507, 

L113 

Yoshii, Y., Tsujimoto, T., & Nomoto, K. 1996, ApJ, 462, 
266 



